Amplitude Noise Supression in Cavity-Driven Oscillations of a Mechanical Resonator 
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We analyze the amplitude and phase noise of limit-cycle oscillations in a mechanical resonator cou- 
pled parametrically to an optical cavity driven above its resonant frequency. At a given temperature 
the limit-cycle oscillations have lower amplitude noise than states of the same average amplitude 
excited by a pure harmonic drive; for sufficiently low thermal noise a sub-Poissonian resonator state 
can be produced. We also calculate the linewidth narrowing that occurs in the limit-cycle states 
and show that while the minimum is set by direct phase diffusion, diffusion due to the optical spring 
effect can dominate if the cavity is not driven exactly at a side-band resonance. 
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Recently there has been considerable progress towards 
the goal of observing signatures of quantum behavior in 
the collective vibrations of mechanical resonators [l[ . Ev- 
idence of quantum behavior is expected to be found in 
the production of non-classical states of mechanical mo- 
tion such as squeezed states or superpositions of spatially 
separated states 0, Q ■ Quantum effects should also be 
evident in the non-linear dynamics of mechanical res- 
onators Q , as well as in the fundamental limits on the 
sensitivity with which mechanical motion can be moni- 
tored 0. However, because mechanical resonators typ- 
ically have resonant frequencies in the radio-frequency 
range or below, thermal fluctuations would naturally 
tend to mask the quantum features and hence significant 
efforts have been devoted to developing ways of cooling 
a mechanical resonator down to its ground state [l|, [H[ . 

The radiation pressure force which arises when a me- 
chanical resonator is coupled parametrically to a driven 
optical cavity [f| provides a very effective way of sup- 
pressing thermal fluctuations in mechanical resonators. 
When the cavity is driven below resonance, quanta are 
absorbed from the mechanical resonator by the cavity. 
The low level of photon noise in the cavity means that, 
provided the relaxation rate of the cavity is much less 
than the mechanical frequency (the good cavity limit), 
the mechanical resonator can in principle be cooled al- 
most all the way to its ground state [6|. However, in 
practice the cooling effect of the cavity competes with 
the resonator's thermal environment and recent experi- 
ments Q have combined the driven cavity with cryogenic 
cooling to achieve lower occupation numbers. 

If the cavity is instead driven above resonance then en- 
ergy is absorbed by the mechanical resonator leading to 
states of self-sustaining oscillation |{J [ljj ■ Increasing the 
power of the cavity drive leads eventually to a region of 
multistability marked by a sequence of dynamical tran- 
sitions between limit-cycle states of different sizes 0, [l(| ■ 
Little attention has so far been devoted to studying the 
quantum aspects of the limit-cycle dynamics, although 
recent numerical calculations began to explore the behav- 
ior in this regime [14(. However, similar laser- like states 



range of finite-level systems [11- 13|. 

In this Letter, we present an analytic calculation of 
the amplitude noise of a cavity-driven mechanical res- 
onator within a limit-cycle. Our principal finding is that 
the amplitude noise in a limit-cycle can be very low: for 
very low thermal noise the resonator can be driven into 
a sub-Poissonian state by the cavity. More generally, at 
a given temperature the amplitude noise in a limit-cycle 
state can be substantially lower than in an equivalent one 
produced by a perfect harmonic drive. We also explore 
the behavior of the resonator linewidth in the limit-cycle 
state, generalizing a previous calculation (l7j . 

The parametrically coupled driven cavity and median 
ical resonator system is described by 
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H = -h[A+ |(6 + fe t )]a t a + %uj m tfb + Ml{a + a)), (1) 

where g is the coupling strength, a(b) is a cavity (res- 
onator) lowering operator, uj m is the mechanical fre- 
quency and VI parameterizes the strength of the laser 
drive. The cavity is driven at a frequency Ud detuned 
from the cavity frequency, w c (3> oj m ), by A = cud — oj c . 
The evolution of the system is described by the master 
equation |14[ , 
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where the coupling of the mechanical resonator to its 
thermalized surroundings at temperature T is described 
by, 
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have been studied in mechanical oscillators coupled to a 



with 7 m the mechanical damping rate and n = 
[exp(hu! m /k-QT) — 1] . The cavity dissipation is de- 
scribed by, C c p = —j c (a^ap + pa^a — 2apa T )/2 with *y c 
the decay rate and we assume Tiuj c k^T so that we can 
neglect thermal fluctuations. 

We proceed by carrying out a Wigner transformation 
of the master equation, which introduces the complex 
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variables a and (3 for the phase space of the cavity and 
resonator respectively [1] . Neglecting third-order deriva- 
tive terms (truncated Wigner function approximation) in 
the resulting equation of motion for the Wigner function 
leads to a standard Fokkcr-Planck equation from which 
we obtain the coupled Langevin equations, 

d = i[A+^(P + p*)]a-iSl-^a + ri a (3) 
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The stochastic force terms [8[ rj a ,r]p have zero means 
and non-zero second order moments (r] a * (t)rj a (t')) = 
<J(t-f )7c/2 and (^.(fjf^f )) = 5(t-t') lm (n+ i). The 
truncated Wigner function approximation is expected to 
describe small linear fluctuations [1] and hence should 
provide a good description for the limit-cycle states. 

We follow the approach used by Marquardt et al. [l^ 
to solve the corresponding classical dynamics and ex- 
tend this to include the noise. We make the realistic 
assumption that the total resonator damping is much 
lower than the cavity decay rate so that the amplitude 
and phase of the resonator change only very slowly on 
the time-scale of the cavity dynamics. The problem 
is now split into two parts. First, we solve for a us- 
ing the ansatz /3(t) = /3 C + Be^e"^™*. The result- 
ing solution is split into average and fluctuating parts, 
a(t) = (a(t)) + 5a (where the average corresponds to 
the solution obtained when the stochastic force term is 
dropped) and then assuming weak fluctuations we ap- 
proximate aa* ~ (eta*) + 8a* (a) + Sa(a*} to obtain an 
effective equation of motion for j3(t) = (3(t) — /3 C . 

Solving for the cavity dynamics and taking the Fourier 
transform, we obtain 



FIG. 1: (Color online) Resonator Fano factor F and average 
energy (n) (inset) as a function of A calculated using P(B) 
(dashed red curves) and numerically (solid blue curves). In 



u) m = 1, g = 0.4, 7„ 



5 x 10" 



and in (b) cj m = 5, 



g — 1.5, 7 m = 3 x 10~ J , in each case n = and Q = 0.05 (we 
adopt units such that j c = 1). The results obtained from the 
Gaussian approximation to P(B) [see text] are shown in (b) 
(dashed green curves). 



which, although non-linear, is well approximated by its 
linear form for weak g. 

We focus for now on the fluctuations in the amplitude 
of the resonator motion. The equation of motion for the 
amplitude can be written as, 
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where -y T = 7m + 1b a and rj^ = i( I | |3 e' ( * +u "' i) ± 
V)3te -i(<p+u> m t)j + if(( a *)5a + (a)5a*)(c l ^ +UJ ^ =p 
e -i(<t>+u m t)y Xhc term rj^ relates to the phase diffusion, 
as discussed below. 

An effective diffusion constant valid on timescales long 
compared to uj m and j c is obtained (ljj from the zero- 
frequency component of the correlator {rj^,(t)rj^,(t')), av- 
eraged over a mechanical period to eliminate explicit time 
^ J2„ J- n {z)e l<l > n 6(LL> - io m dependence, 



7c/2 + i(uj - A) 



(5) 

and 5a'{ui) = i]' a /[y c /2 + i(uj — A)], where z = gB/uj m , 
J n (z) is a Bessel function of the first kind, A = 
A + gRe[/3 c ] and the primes denote e.g. rj' a (t) = 
ri a (t)e~ lz Sla ^ +wt ' . Keeping only the fundamental oscil- 
lating component of (a) (a* ) , we obtain 

(6) 

where the effective damping and frequency shift of the 
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The contribution from the resonator's thermalized sur- 
roundings is J D m = \"f m {n + |), and the contribution 
from the cavity is given by, 
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resonator due to the cavity are given by [10( , 
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with h r . 



■ + i(A + nuj m ). The center of the mechanical 



oscillations is given by, 



Dg A again relates to the phase diffusion. 

The Fokker-Planck equation equivalent to Eq. ((9|) has 
a steady-state solution P(B) oc cxp(— U(B)), with 
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where small corrections to the drift terms due to the noise 
have been neglected [IB, 16 1. This potential solution can 
be used to calculate the average resonator energy, (n) 
(n = b'b), and associated fluctuations over a wide range 
of parameters, including both the good and bad cavity 
limits. Figure [T] shows a comparison of (n), and the res- 
onator Fano factor F = ((n 2 ) — (n) ) / (n), obtained using 
the P(B) distribution and the results of a direct numer- 
ical solution [ij of the master equation [Eq. ([2])]. The 
numerical calculation is performed in a restricted num- 
ber state basis using 3 states for the cavity and up to 120 
for the resonator. We also neglect elements representing 
coherence between resonator states with a large separa- 
tion in energy [l2j]. Representing the cavity in a basis 
centered on the equivalent uncoupled (g — > 0) state, i.e. 
a' = a + iO/(7 c /2 — iA), allowed us to study strongly 
driven cavities using only a few states (so long as the 
associated variance is not too large). 

We see from Fig. [TJ that there is very good agreement 
between the numerics and the calculation using Eq. ([12")) 
when the resonator is in a limit-cycle state (character- 
ized in Fig. Q]by a large value of (n) and a relatively low 
Fano factor.) The agreement is still quite good [l9| when 
the resonator undergoes a dynamical transition from the 
limit-cycle state back to one in which it fluctuates instead 
about a fixed point [ll|, El, EH (marked by peaks in the 
Fano factor). The main approximation we have made is 
to neglect the higher-order derivatives (and hence higher- 
order correlations) by truncating the Wigner function. 
The relatively strong couplings we used in order to cap- 
ture the limit-cycle dynamics numerically [bj ] provide a 
severe test of this approximation. The slight shift be- 
tween the analytical and numerical curves in the figures 
is a sign that we are approaching the limits (in terms of 
coupling strengths) of the validity of this approach [li| . 

In a limit cycle the resonator distribution can be ap- 
proximated as a Gaussian centered at an amplitude Bq, 
determined by the condition 7 m = — jba{Bo), with a 
width given by a 2 = (D rn + D BA (B ))/(2B ^\ Bo ), 
and a Fano factor F ~ 4a 2 . The Gaussian approxima- 
tions to (n) (given by B 2 — 1/2) and the Fano factor are 
compared with results from numerics and using the full 
P(B) distribution in Fig. [TJ). 

The resonator Fano factor drops when it is in a well- 
defined limit-cycle state. In the highly idealised case 
where thermal noise is negligible (n s=s 0), the cavity can 
drive the resonator into a non-classical sub-Poissonian 
state with F < 1 (sec Fig. [2]). This is a consequence of 
the low noise properties of the cavity: it is well-known in 
the context of laser physics that regular pumping can 
lead to sub-Poissonian states 0]. However, for a me- 
chanical resonator thermal noise plays an important role 
and we now examine to what extent we can think of 
this being suppressed in the limit-cycle states. A useful 
comparison can be made between the cavity-driven res- 
onator states and a displaced thermal state (DTS) with 
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FIG. 2: (Color online) Analytic [Eq. (fl2|) ] (dashed red curves) 
and numerical (solid blue curves) calculations of the resonator 
Fano factor F in terms of Fd- The curves are (top to bottom) 
n = 0,0.5, 1,5. In each case w m = 5, r y m = 3 x 10~ 5 , g = 2 
and Q = 0.05. The dots represent the on-resonance value 
calculated using Eq. (|13|) . Note, for n = 0, Fd — 1 and F 
reaches a minimum ~ 0.9. 



the same amplitude, thermal occupation number n and 
external damping j m . A DTS [2(| is produced by har- 
monically driving a resonator initially in a thermal state, 
increasing its energy without introducing additional fluc- 
tuations. We choose a DTS for comparison as it reduces 
to a coherent state for n — > 0, meaning it is both the gen- 
eralization of a coherent state to finite n and of a thermal 
state to finite amplitude. 

The Fano factor of a DTS with amplitude Bq is Fd = 
(n(l + n) + (2n + 1)\B \ 2 ) / (fi + \B \ 2 ) and we compare 
this with that of the cavity-driven resonator in Fig. [2] for 
a range of external bath temperatures. In a well-defined 
limit cycle the value of the ratio F/Fd is suppressed sig- 
nificantly below unity and we can think of the cavity as 
suppressing the thermal fluctuations. Note that this sup- 
pression also occurs outside the good cavity limit shown. 
As with the usual cavity induced cooling the noise 
suppression occurs because the low-noise cavity can in- 
crease the friction (damping) on the resonator without 
adding significantly to the diffusion. Thus, driving a me- 
chanical resonator via a cavity in this way produces a 
lower-noise state than driving with a perfect harmonic 
drive. 

For A ~ u) m and in the good cavity limit, we can use 
Eqs. and ([T2"]). together with the fact that P(B) is 
almost Gaussian in a well-defined limit-cycle to obtain a 
simple approximate expression for the Fano factor, 
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where the amplitude Bq is defined by "f rn = —^iba{Bq)] 
the predictions of this equation are shown as dots in Fig. 
[2] This formula breaks down when the cavity-resonator 
coupling is increased sufficiently to allow the co-existence 
of more than one stable limit cycle [l(| , (although it cor- 
rectly predicts F within the second limit cycle once the 
first has become unstable). Thus whilst Eq. (fTB")) sug- 
gests that an arbitrarily small F can always be achieved, 
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the actual minimum value achievable for a given system 
is set by this expression together with the requirement 
that 7 m = —^ba(B) has a single (non-zero) solution. 

As well as fluctuations in amplitude, the resonator also 
undergoes phase diffusion which determines the linewidth 
in the limit-cycle state. In a well-defined limit-cycle, we 
can write a coarse-grained equation of motion for the 
phase 11, lH </>, again linearizing the fluctuations, 



Soj l SB. 



(14) 



where SB represents the amplitude fluctuations B — Bq 
and Sujl = Bo^§-\b is the frequency shift linearized 
about the limit cycle. Defining the phase diffusion in the 
same way as the amplitude diffusion, Eq. (|10[) , we get, 
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\b is the linearized damping and 



cross-correlations have been neglected. Because the shift 
in the resonator frequency due to the cavity (optical 
spring effect) is amplitude dependent, the amplitude fluc- 
tuations can give rise to an important additional con- 
tribution to the phase diffusion. The phase diffusion is 
shown in Fig. [3] and it is clear that although the opti- 
cal spring contribution is negligible at the center of the 
side band resonance (where 8u>l itself is negligible 17j). 
it nevertheless becomes important on either side. 

Numerical calculations of the resonator spectrum 
S(uj) = J dte- luJt {{tf(t),b(0)}) show ncar-Lorentzian 
peaks around uj = and lj = u) m of width Ao and A Wm . 
These widths are determined by the slowest dissipativc 
timescale at each frequency [l2 , 15 1 , and hence Ao is given 
by the energy relaxation rate A = ^ 1% A I b anc ^ ^he 
linewidth is given by the phase diffusion, K UJm = 
There is good agreement between the numerical and an- 
alytic calculations of these quantities (Figs. |3]b,c) within 
the limit cycle regime (the main difference again an ef- 
fective shift in A). 

In conclusion, we have studied the amplitude and 
phase noise of limit-cycle states of a mechanical res- 
onator driven by an optical cavity. Within a limit-cycle 
amplitude fluctuations are suppressed in the sense that 
they can be substantially less than in a corresponding 
state produced by simply applying a pure harmonic drive, 
the counterpart of the cooling that occurs in the stable 
regime. For low enough thermal noise the cavity gener- 
ates non-classical sub-Poissonian resonator states. How- 
ever, for the phase diffusion in the limit-cycle states the 
cavity noise simply adds to the effects of thermal fluctu- 
ations and the optical spring effect can also generate a 
significant contribution. The quantum noise of the res- 
onator is described rather well by the truncated Wigner 
function approach over a range of resonator frequencies, 
both within the limit-cycle states and more surprisingly 
within the transition regions. 
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FIG. 3: (Color online) a) Phase diffusion calculated nu- 
merically (solid blue) and analytically (dashed red). Solid 
and dashed thin lines show the direct phase diffusion 
(Dm + Dg A )/2 and that due to the optical spring effect 
28u\/'y\(D rri + Dg A ) respectively. The insets show the 
linewidths Ao (b) and A Um (c). Parameters are: u m = 
5,7 m = 3 x 10" 5 ,# = 1.5, fi = 0.05, n = 0. 
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